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Abstract 

We study the structure and reveal the physical nature of the reversed granulation pattern 
in the solar photosphere by means of 3-dimensional radiative hydrodynamics simulations. We 
used the MURaM code to obtain a realistic model of the near-surface layers of the convection 
zone and the photosphere. The pattern of horizontal temperature fluctuations at the base of 
the photosphere consists of relatively hot granular cells bounded by the cooler intergranular 
downflow network. With increasing height in the photosphere, the amplitude of the tempera- 
ture fluctuations diminishes. At a height of z — 130 — 140 km in the photosphere, the pattern 
of horizontal temperature fluctuations reverses so that granular regions become relatively cool 
compared to the intergranular network. Detailed analysis of the trajectories of fluid elements 
through the photosphere reveal that the motion of the fluid is non-adiabatic, owing to strong 
radiative cooling when approaching the surface of optical depth unity followed by reheating 
by the radiation field from below. The temperature structure of the photosphere results from 
the competition between expansion of rising fluid elements and radiative heating. The former 
acts to lower the temperature of the fluid whereas the latter acts to increase it towards the 
radiative equilibrium temperature with a net entropy gain. After the fluid overturns and de- 
scends towards the convection zone, radiative energy loss again decreases the entropy of the 
fluid. Radiative heating and cooling of fluid elements that penetrate into the photosphere and 
overturn do not occur in equal amounts. The imbalance in the cumulative heating and cooling 
of these fluid elements is responsible for the reversal of temperature fluctuations with respect 
to height in the photosphere. 

* Present address: Lockheed Martin Solar and Astrophysics Laboratory, Bldg/252 3251 Hanover St., Palo Alto, 
CA 94304, USA. 



1 Introduction 



In the quiet Sun, the most prominent photospheric feature is the surface granulation pattern. 
Observations in white light show bright granules with a typical size of 1 Mm, embedded in a network 
of dark intergranular downflow lanes and vertices. The granulation pattern is non-stationary: 
studies of time-sequences of the granulation pattern in t he quiet Sun show that individual features 



evolve with an average lifetime of about 5 minutes (e.g. Title et al. . 19891 ). 



While intensity images in the visible continuum taken near disk-center show bright granules 
and comparatively dark intergranular boundaries, observations taken in the wings of the Ca II H 
fc K lines exhibit a similar pattern, but wit h reversed intensity contrast (jEvans fc Catalano . 
1972 ; Suemoto. Hiei. fc Nakagoml 1987 . 19901 ). This effect, known as the reve rsed granulation , 
has also been detected in intensity measurements of photospheric spectral lines (lEspaenet et al 



1995; Balthasar et al. , 1990; iKucera et alj . 1995; iRodriguez Hidalgo et al.l . [l999l ; |Janssen fc Cauzzil 



2006). These observations suggest that the reversal of the intensity contrast occurs already at a 
height of abo ut 150 km above the mean geometric level where T500 = 1. 

Recently, iRutten. de Wijn. fe Sutterhnl (120041 ) acquired simultaneous, speckle-reconstructed 
images of quiet-Sun granulation in the G-band and the Ca II H line. On the basis of cross- 
correlations between time sequences of intensity images in these two spectral ranges they found 
that the anticorrelation between the two patterns is highest if the Ca II H intensity maps are time- 
delayed by 2 — 3 minutes with respect to the G-band intensity maps. By computing anticorrelation 
coefficients between (1) intensity images in the continuum and in the Fe I 7090 A line and (2) i nten- 
sity images in the continuum and in the wing of the Ca II 8542 A line, Janssen fc Cauzzi ( 20061) also 
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found that the strongest antic orrelation occurs for a time delay of ~ 2 min. 
(|2005[ ) extended the work of IRutten. de Wiin. fc Siitterlin ( 20041 ) by calculating synthetic inten- 
sity images in the blue continuum and in the Ca II H wing (in the LTE approximation) from a 
3-dimensional radiative hydrodynamic simulation. From the unsmoothed synthetic intensity im- 
ages, they found that the anticorrelation between the continuum images and the Ca II H wing 
images does not vary significantly between a time delay of and 2 minutes. The anticorrelation 
decreased for longer time delays. When the images were smoothed to 1.5 arcsec resolution, a time- 
delay of 2 — 3 minutes yielded a significant l y high er value of the anticorrelation compared to zero 
time delay. IRutten. de Wiin. fc Sutterhnl (2004) speculated tha t magnetic fields play no major 
role in the formation of the reversed granulation. The results of Leenaarts fc Wedemever-Bohml 
(2005) is evidence supporting this view. 

Numerical simulations have firm l y established the solar granulation as a radia tive-convective 
phenomenon jNordlundl . 1 1984 1 19851: isteffen et all 1 19891: 1 Stein fc Nordlundl . Il998j) . Granules are 
upwellings of hot, high-entropy material originating from the convection zone and overshooting into 
the stably stratified photosphere. Radiative losses make the overshooting material relatively cold 
and dense. Its overturning motion supplies the network of intergranular downflows. Simulations 
have also revealed that the pattern of horizontal temperature fluctuations in the middle/upper 
photosphere is inverted compared to the corresponding pattern at the base of the photosphere. 
Whereas at the base of the photosphere, the granules are relatively hot compared to the intergran- 
ular network, the former are relatively cooler for T500 < 0.3. It is well known that this reversal in the 
temperature fluctuations in the photosphere is responsible for the reversal in the observed inten- 
sity p attern ( Nordlundl . 19841 : ISuemoto. Hiei. fc Nakagoml 1987 : Leenaarts fc Wedemever-Bohm . 
20051) . 



What is less clear is the physical mechanism which leads to the reversal of the temperature fluc- 
tuations with height. It is sometimes incorrectly assumed that flui d travelling in th e optically thin 
layers of the photosphere does so adiabatically. As pointed out bv lNordlundl (|1984l ). this cannot be 
true because material with an original temperature of 5900 K penetrating and expanding adiabat- 
ically in the photosphere would reach a temperature of 1800 K after traversing only three pressure 
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scale heights. Such low temperatures are not detected in the photosphere. IStein Sz. Nordlund 
(|1998T ) indeed find that fluid rising in the photosphere is typically reheated by the radiation field. 
This leads one to conclude that non-adiabatic effects, such as radiative energy exchange, must play 
a crucial role in the formation of the reversed temperature fluctuations. 

The aim of the present paper is to investigate the physical origin of the reversal of temperature 
fluctuations. As we shall see in the following sections, radiative energy exchange is indeed key to 
the temperature fluctuation reversal. By following in detail the time histories of fluid elements 
penetrating into the photosphere, we identify how radiative heating and cooling, in tandem with 
the horizontal pressure fluctuations, causes the reversed granulation. The paper is structured as 
follows. In § [21 we describe the set of governing equations, the numerical methods used and the 
simulation setup. In § [3j we discuss results from the numerical simulation, including the mean 
stratification of the 3D model (§ 13. ip . the temperature structure in the photosphere (§ 13. 2[) and 
reveal the physical origin of the reversed granulation (§ 13. 3p . Finally, in § HI we discuss how the 
present work dispels previous misconceptions and improves our understanding of the phenomenon. 



2 Simulation setup 



2.1 Governing equations 

For a realistic treatment of the dynamics and energetics in the near-surface layers of the non- 
magnetic solar atmosphere, it is important to co nsider the effects of compressibility, energy ex- 
change via r adiat i ve trans f er and partial ionization ([Stein fc Nordluridl .fi 998: Ra st fc ToomreL 1993 : 
iRast et aL . 19931 : Vogler . 2003 ). The hydrodynamics equations, the radiative transfer equation 
(RTE) and the equation of state (EOS) including the effects of partial ionization constitute the 
system of governing equations for our model. Since we are confining our attention to the near- 
surface layers of the convection zone and the photosphere, the approximation of local thermody- 
namic equilibrium (LTE) suffices for the treatment of the radiative transfer. The hydrodynamics 
equations include the continuity equation, the momentum equation and the energy equation, viz. 
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(3) 



where v is the fluid velocity, g the mass density and e the internal energy density per unit mass. 
The operator (g) denotes the tensor product. In the near-surface layers of the Sun, the gravitational 
acceleration is almost constant and we take g = —2.74 x 10 4 z cm s~ 2 . 

The hydrodynamics equations ([I])-© must be supplemented by the appropriate constitutive 
relations describing the material properties of the gas. These include the equation of state (EOS), 
functions specifying the viscous stress tensor r and thermal conductivity K, as well as functions 
specifying the radiative properties of the gas (e.g. source function S u , opacity k„ etc.). The 
frequency-integrated radiative heating term appearing on the r.h.s. of Eq. |3]) is given by 



Qr 



ad 



(4) 



where J„ is the angle-averaged specific intensity and S v the angular-averaged source function. 
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2.2 Numerical method 



For our radiative hydrodynamics simulations, we have used the MPS/ C/ni versify of Chicag o 
Radiative MHD (MURaM) code (|Vogler. Shelvag. Schiissler. Cattaneo. Emonet. fc Lindel . 120051 ). 
In the absence of magnetic fields, the equations solved by MURaM reduce to the equations given 
in the preceding section. 

In the following, we briefly desc ribe the nume rical methods implemen ted in MURaM. For 
further details, we refer the reader to Vogler ( 20031 ) and Vogler et al.l ( 20051 ). MURaM integrates 
the discretized equations ([I])-® on a 3-dimensional cartesian grid, with constant grid-spacing 
in each cartesian direction. The spatial discretization of the equations is based on a fourth- 
order centered-difference scheme on a five-point stencil. Time-stepping is carried out using an 
explicit fourth-order Runge-Kutta scheme, with the maximum size of the time step determined by 
a modified CFL-criterion. 

In order to prevent the nonlinear cascade of energy towards the grid scale, artificial diffusivities 
are implemented to stabilize the s imula tions. A combination of hyper- and shock-resolving diffu- 
sivities are used (jCaunt fc Korpil . l200lh . Hyperdiffusivities are implemented in such a way, that 
the diffusion terms associated with them are comparable in magnitude to the inertial terms only 
at length scales close to the grid-spacing. For the thermal and viscous c oeffici ents K and v, the 
hyperdiffusivities are defined using the scheme described in IV5gler et al. I (|2005l) . 

The RTE is solved with the short-characteristics method on a 3D cartesian grid, with the 
grid centers of this radiative grid defined by the cell corners of the grid used for discretizing the 
hydrodynamics equations ([I])-©- For each point on this radiative grid, the RTE is solved over 24 
ray directions (3 per octant). The first moment of I v over these rays provides the radiative flux 
density F^, which in turn allows Q ra d to be calculated. 

Voglerl liool and 



see 



MURaM can treat the radiative transfer in the grey or non-grey case 
references therein). For this study, we have carried out two simulations, namely a grey run and 
a non-grey run. The results presented in the rest of this paper correspond to the non-grey case. 
The main difference between the results of the two cases it that the non-grey case yields smaller 
temperature fluctuations in the photosphere than the grey case. This difference arises because the 
grey case underestimates the smoothing effects of radiati ve tran s fer in the optically thin layers of 
the photosphere in comparison with the non-grey case (jVoglerl . |2004) . Our explanation for the 
reversed granulation phenomenon (see § 13. 2p . however, applies to both cases. 

The EOS relations T — T(e, g) and p = p(e, g) are used by MURaM in the form of look-up 
tables. The look-up tables were evaluated by solving the Saha-Boltzmann equations. The first 
ionization of the 11 most abundant elements was considered. For a detailed account of how the 



tables were evaluated, we refer the reader to Appendix A of IVogler et alJ (|2005l ) 



2.3 Simulation setup 

In order to obtain a thermally relaxed 3D model of the upper convection zone and the photosphere, 
we have chosen the following setup. The horizontal size of the simulation domain is 24 Mm xI2 
Mm, with each direction spanned by 480 and 240 grid cells, respectively. The height of the domain 
is 2.3 Mm, and is spanned by 144 grid cells. The grid-spacing in the horizontal and vertical 
directions is 50 and 16 km, respectively. Periodic boundary conditions are assumed for the vertical 
side boundaries. The bottom bound ary is open and allows mass to flow smoothly into and out of 
the domain (see Vogler et al. . 20051 ). 

For t he initi a l con dition, we took the ID height profiles of g and e from the mixing length 



model of Spruit! ( 1974 ) and introduced a plane-parallel stratification into the simulation domain. 
To break the symmetry, we imposed pressure perturbations on the order of a few percent. We then 
ran the simulation until we obtained a thermally relaxed, statistically stationary state. 
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Figure 1 : Horizontally averaged pressure (solid) and temperature (dashed) as functions of height 
z. 

3 Simulation results 

3.1 Mean stratification 

The horizontally averaged profiles of pressure and temperature as functions of height are shown 
in Fig. [TJ respectively, as solid and dashed curves. The geometrical height z = corresponds 
to the geometrical level where the continuum optical depth at 500 nm is, on average, unity. 
T he height dependence of t he mean profiles are in general agreement with the simulation results 
of Stein & Nordlundl (Il998h . 



Figure [21 shows the logarithmic temperature gradients determined from the mean profiles: 

V := — , and 5 

dmp 



' ad 



91nT 
dlnp 



(6) 



V (dotted curve) is the actual average temperature gradient in the simulation and V a d (solid black 
curve) is the adiabatic temperature gradient. V a d describes the variation of temperature of a fluid 
element undergoing adiabatic expansion or compression. 

The stability of the stratification to convective motion is affected by changes in the ionization 
state. The superadiabaticity of the stratification is defined as S := V — V a d and is plotted as a 
dashed curve in Fig. [2] In the photosphere, the chemical species are almost completely monatomic. 
For such a gas, 7 = 5/3 and hence, V a d = 1 — 1/7 = 0-4- With increasing depth into the convection 
zone, the ionization fraction of hydrogen (and traces of other species) increases. This has the effect 
of decreasing the adiabatic temperature gradient, so that V a( j < 0.4 while 5 remains positive. In 
other words, changes in the ionization state of the gas reduces V a( j and thus enhances buoyancy 
driving. 
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Figure 2: Logarithmic temperature gradients V (dotted curve), V a d (solid black) and the supera- 
diabaticity S = V — V a( j (dashed curve). 



3.2 The reversal of temperature fluctuations in the middle photosphere 

In this section, we study the reversed granulation in our 3D model. For lines that form under 
approximate LTE conditions, such as some Fe I lines and the wings of the Ca II H & K lines, 
fluctuations in line intensity are a proxy for diagnosing temperature fluctuations. The intensity 
image therefore maps the temperature distribution in the layer of the atmosphere in which the 
corresponding part of the line is formed. An inspection of the temperature fluctuations at different 
heights of the atmosphere allows one to determine the height at which the temperature fluctuations 
start to reverse. 

The upper row of Fig. [3] shows the relative temperature fluctuations (AT/T) at the surfaces 
T 500 — 1 (left column) and T500 = 0.1 (right column). The lower row gives the vertical velocity at 
the same surfaces. On average, the T500 = 0.1 surface is elevated above the T500 = 1 surface by 
160 km. At T500 = 1, one finds the normal granulation pattern: the granules are hotter than the 
intergranular lanes. At T500 = 0.1, the reversed granulation is already clearly discernible: the lanes 
are regions of temperature excess. When we inspect the velocity structure, we find that the pattern 
at T500 = 0.1 is very similar to that at T500 = 1, albeit with smaller flow speeds. The upflows in 
the higher layers result from the overshooting of upwelling gas from the convection zone. Thus, 
the reversal of the pattern of temperat ure (and intensity) fluc tuations is not accompanied by a 
reversal in the velocity pattern (e.g. see Ruiz Cobo et al. . ll996l) . 



In order to determine the geometrical height at which the reversed granulation begins, we 
consider the vertical velocity and temperature as functions of z for vertical lines of sight (l.o.s.). 
Distinguishing the l.o.s. profiles according to the sign of the vertical velocity at T500 = 1, we 
separately determine averages for profiles of the temperature fluctuations in granules and in the 
intergranular downflow network, respectively. The fluctuations are evaluated relative to the overall 
average temperature profile. These two quantities are plotted as functions of z in Fig.|¥J Diamonds 
indicate the mean temperature fluctuation in the granules whereas crosses indicate the mean tem- 
perature fluctuation in the intergranular lanes. These two are not simply the negative of each 
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Figure 3: Temperature fluctuations (AT/T, upper row) and vertical velocity (lower row) at two 
surfaces of constant optical depth. The left column shows the patterns at T500 = 1, the right 
t~5qo = 0.1. The reversed granulation corresponds to a reversal of the temperature fluctuations 
while the velocity pattern remains qualitatively unchanged. 
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Figure 4: Mean temperature fluctuation (AT) of granules (diamonds) and of intergranular lanes 
(crosses) as a functions of geometrical height. At z = 130 — 140 km, the cells and boundaries 
have, on average, the same temperature. Above this geometrical height, the intergranular lanes 
are relatively hotter than the granules. 



other because the cells have a larger area fraction than the boundaries. In the convection zone 
and in the lower photosphere, upflowing material in the granules are relatively hot compared to 
downflowing material in the lanes. The reversal of the horizontal temperature fluctuations occurs, 
on average, in the middle photosphere at a height of z = 130 — 140 km. Above that height, the 
intergranular lanes are hotter than the granules with a temperature excess of 100-200 K. 

We have also determined the temperature profiles in granules and in the intergranular network 
as functions of optical depth along vertical lines of sight. The average temperature fluctuations 
of the granules and the intergranular lanes as functions of T500 are indicated by diamonds and 
crosses respectively in Fig. [5] In the convection zone (z < 0), the temperature fluctuations at a 
given geometrical height are much larger than th e temperature fluctuations at constant optical 
depth. As pointed out bv lStein fc Nordlundl (l998), this is due to the high temperature sensitivity 



of the H~ contribution to the continuum opacity. In this figure, the two curves cross each other 
a t T500 = 0.22. On average, the surface of constant optical depth at T500 = 0.22 is elevated above 
the T500 = 1 surface by 120 km. 

The fact that the two different averaging procedures yield somewhat different heights (z = 
130 — 140 km and z — 120 km, respectively) for the layer of the photosphere where the tem- 
perature fluctuations reverse their sign is easily understood. In the former case, we calculated 
the horizontal temperature fluctuations relative to the mean temperature of the layer at differ- 
ent geometrical heights. In the latter case, the temperature fluctuations correspond to deviations 
from the mean temperature at surfaces of constant optical depth. Since these surfaces are cor- 
rugated (e.g. the rms fluctuation of the geometrical height of the T500 = 1 surface is 30 km), 
temperature fluctuations across such a surface only approximate the horizontal temperature fluc- 
tuations at the same geometrical height. In terms of both geometrical height and optical depth, 
our determination of the level at which the reversal occurs is in good agreement with previous ob- 
servational work ( Ruiz Cobo et all Il996& iRodriguez Hidalgo et all Il999l : iBorrero k. Bellot Rubiol . 
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Figure 5: Mean temperature fluctuation (AT) of granules (diamonds) and of intcrgranular lanes 
(crosses) as a functions of optical depth. The reversed granulation appears at optical depths of 
about T500 = 0.2. 



2002t IPuschmann et all l2005h . 



3.3 The origin of the reversal of temperature fluctuations 

The analysis in the previous section confirmed the reversal of horizontal temperature fluctuations 
in the middle photosphere. Now we proceed to explain the physics behind this phenomenon. 
There are two important aspects to the explanation. Firstly, we need to understand the nature of 
convective overshoot of granular fluid into the photosphere: namely, what is the balance of forces 
controlling the dynamics of convective penetration and overturn? Secondly, we need to understand 
the thermodynamical processes experienced by the overshooting material. 

In order to address these issues, we have examined the histories of tracer fluid elements in 
our simulation. At time t — (corresponding to the snapshot shown in Fig. [3]), tracer fluid 
elements were introduced just beneath the surface inside a granule and propagated according to 
the simulated velocity field. The physical properties of the fluid elements (e.g. v, g, e, s, Q ra d etc.) 
at th eir actual positions were d etermined by tricubic interpolation from the values in the adjacent 
cells (Le kien fc Marsdenl . [2005h . 



Figure \E\ shows a subvolume of the simulation domain in the neighbourhood of the granule 
centered at [x, y] — [7, 8] Mm in Fig. [3J The greyscale image at z = indicates the vertical velocity 
at z = at the instant when the tracer fluid elements were released (t = 0). The colour-coding 
of the trajectories indicate the sign of Q ra d- dark blue means Q ra d < (cooling), yellow means 
Qrad > (heating). The set of trajectories resembles a fountain: the fluid elements originate from 
the granule interior, travel up to the photosphere and eventually overturn and descend into the 
downflow network in the convection zone. 
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Figure 6: Trajectories of tracer fluid elements in the photospheric neighbourhood of a granule. The 
greyscale image in the z = plane shows the vertical velocity there at the time when the tracer 
fluid elements were released. The colour coding of the trajectories corresponds to the sign of the 
radiative heating Q ra d experienced by the fluid elements (dark blue means Q ra d < 0, yellow means 
Qrad > 0). 



3.3.1 The motion of fluid elements in the photosphere 

First we concentrate on the motion and dynamics of the fluid elements. As Fig. [5] shows, in 
spite of the strongly stable average stratification, the fluid elements are able to penetrate into 
the photosphere. The shape of the trajectories tell us that in the photosphere, they experience a 
deceleration in the vertical direction and a deflection towards the edge of the granule. 

In the following, we give a brief description of the force balance which drives the horizontal 
granular outflows and the feeding of the downflows in the intergranular lanes. Inspection of the 
horizontal pressure fluctuations in the photospheric layers of our simulation reveal that local pres- 
sure hills in the interior of granules are responsible for the buoyancy braking of the upflows and the 
deflection of material horizontally outwards from the centres of granules. At the interface between 
two granules, the horizontal outflows from the granules meet and develop localized pressure ridges. 
These localized regions of pressure enhancement serve to retard the horizontal granular outflows. 
So, when a fluid element in the granular outflow approaches the intergranular boundary, it is decel- 
erated by the pressure gradient. Finally, the pressure enhancement at the boundaries, together with 
gravity, work in tandem to accelerate the fluid element downwards, thereby feeding the intergran- 
ular downflow network. Th ese results are in accordance with results from 2D idealized simula tions 
of compressible convection ( Hurlburt k. Toomre . 19881: iHurlburt. Toomre. fc Massaguer , 1986f ) and 



3D near-surface convection simulations ( NordlundT 19851 : IStein fc Nordlundl . 1998| ). 



3.3.2 The thermodynamic history of fluid elements in the photosphere 

To avoid confusion in the following discussion, at this point we would like to clarify the following 
terminology: By heating (cooling), we refer to a non-adiabatic process which increases (decreases) 
the specific entropy of a fluid element. For example, when a fluid element absorbs more radiation 
than it emits, it is heated and gains entropy. A fluid clement which is compressed adiabatically, 
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however, will have a higher final temperature but is not heated. In the absence of viscous dissipation 
or thermal conduction, the equation governing the evolution of the specific entropy of a fluid element 
is 

DS _ Qrad 

Dt ~ qT ' 



(7) 



where D/Dt denotes the Lagrangian derivative (jMihalas fe Mihalad . fl985h . The corresponding 
equation governing the temperature of a fluid element is 



DlnT 
Dt 



(7: 
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Ding 



Dt 

( 73 - 1)V- v + 



dlnT 

ds 
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Ds 
Dt 

Qrad 



(8) 
(9) 



where 73 = 1 + ( gi"^ > 1 is Chandrasekhar's third adiabatic exponent. Two effects contribute 

to changes in the gas temperature of a fluid element. The first term on the r.h.s. of Eq. (Jl]) 
describes the role of adiabatic expansion or compression. Clearly, the effect of adiabatic expansion 
(V • v > 0, Q ra( j = 0) is to decrease the gas temperature. The second term on the r.h.s. of Eq. J9J 
describes the role of radiative heating. As expected, radiative heating (cooling) leads to an increase 
(decrease) of the gas temperature. The relative importance of the two effects in controlling the 
fluid temperature depends on their characteristic timescales. Inspection of Eq. © motivates us 
to define the dynamical and radiative timescales respectively as 



^dyn 
^rad 



1 



(73 - 1)V • v 
ds \ gT 



and 



dlnT J g Q rad 



(10) 
(11) 



In the case that idyn *C i r ad, the motion of a fluid clement is essentially adiabatic. In the op- 
posite extreme, where idyn iS> ^rad> the fluid element is always in radiative equilibrium with its 
surroundings. 

Now we proceed to discuss the thermodynamic history of fluid elements overshooting into the 
photosphere. Owing to the large opacity for z < 0, fluid in the upwellings in the convection 
zone ascend almost adiabatically. The decrease of a fluid element's temperature during its rise 
in the convection zone is due to adiabatic expansion. As it approaches the photosphere, the 
opacity drops rapidly and the fluid element loses entropy by radiative cooling. Thereafter, the 
fluid clement overshoots into the stably stratified photosphere. Although the fluid element is now 
optically thin, it is not completely transparent to radiation and continues to exchange energy with 
its surroundings. In regions above granules, (consisting of overshooting material), radiative energy 
exchange predominantly leads to heating. This effect can be seen in the colour-coding of the 
trajectories in Fig. [6l This figure clearly indicates that as the fluid elements penetrate into the 
photosphere (and as they are diverted toward the intergranular lanes), they experience radiative 
he ating. 



Stein fc Nordlundl (|1998I ) also found that rising fluid in the photosphere is systematically re- 
heated by the radiation field (see Fig. 18 of their paper). The reason for this behaviour is the 
thermostatic action of the radiation field to keep the fluid elements close to radiative equilibrium. 
The radiative equilibrium temperature Tre of a fluid clement embedded in a radiation field is de- 
fined as the temperature at which the fluid element emits and absorbs radiation in equal amounts 
(i.e. Qrad = 0). One can make a simple estimate for the lower limit of Tre by assuming the fluid 
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element to be located at r = 0. In the approximation of grey radiative transfer, k v is replaced by 
the mean opacity R, so that Eq. ^ becomes 

Q Iad = 4<KiiQ[J-B(T)}, (12) 

where B(T) = f B v (T)dv = (cr/7r)T 4 is the frequency-integrated Planck function and J = 
(47r) _1 J 47r J °° I^di/dfi is the angle-averaged, frequency-integrated intensity. Assume that the solar 
atmosphere has plane-parallel symmetry In this case, l v is symmetric about the z-axis. Making 
use of the Eddington-Barbier approximation, we obtain 

Air 



4tt 7 
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(T| T=AI ) 4 d M , (13) 
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where fj, = \± ■ z = cos 9. Here /i is the direction cosine vector and 9 is the angle between \x and 
z. By definition, Tre is such that -B(Tre) = ( cr / 7r )T , RE = J- To evaluate J from Eq. (TTBl . 
we use the average temperature profile T(t) from our model. This yields Tre = 4600 K. This 
result tells us why granular material ascending in the photosphere must be radiatively heated. 
Consider a fluid element that is rising in the photosphere. If the fluid element rises adiabatically, 
its temperature would drop so quickly with height, that it becomes cooler than the radiative 
equilibrium temperature. As soon as this happens, the fluid element will absorb more radiation 
than it emits (Q ra d > 0). By the same token, consider a fluid element in the photosphere which 
has been transported from the interior of a granule towards an intergranular lane. As it descends, 
it is compressed. If it were to descend adiabatically, its temperature would soon be above Ire, 
which means it must then cool radiatively. 

Figure [7] provides an illustration of how overshooting fluid elements become cooler than Tre as 
they penetrate into the photosphere, and subsequently become hotter than Tre as they descend 
towards the convection zone. In the figure, the height z of a fluid element is plotted as a function 
of the horizontal displacement from its original position (i.e. [Aa; 2 + Ay 2 ] 1 / 2 ). In both panels, the 
fluid elements travel from the left to the right of the plot. In the upper panel, the trajectories are 
colour-coded according to the sign of Q ra d- In the lower panel, they are colour-coded according to 
the instantaneous, local value of T — Tre- Tre is evaluated by keeping track of Jj and rt, (i refers 
to the index of the frequency bin) along the trajectories and then finding a temperature T such 
that J2 «»B»(T) = J2 K iJi- 

As a rising fluid element approaches the base of the photosphere, its temperature can be up 
to 400 K higher than Tre- This causes the fluid element to radiate intensely. Once the fluid 
element has reached the optically thin layers of the photosphere, adiabatic expansion and radiative 
heating compete to control the temperature of the fluid element. The relative importance of these 
two effects (measured by the ratio of tdyn and t ra d) determine whether the fluid element will be 
driven towards or away from radiative equilibrium. For a fluid element penetrating into the middle 
photosphere in a granular environment, we find that 0.1 < idynArad ^ 1- This means that the 
decrease in fluid temperature due to expansion of the fluid clement is only partially compensated 
by radiative heating. Consequently, its temperature decreases and it departs further away from 
Tre- In Fig. [71 we can see a general trend for T — Tre to become more negative as a fluid element 
rises. For fluid elements that penetrate beyond z = 250 km, T — Tre can reach a value of —200 K. 

Figure [S] shows the profiles of s, T and g along one of the trajectories shown in Fig. [Sj The 
colour-coding indicates the sign of <5 ra d experienced by the fluid element along its trajectory (dark 
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Figure 7: Upper panel: Trajectories are colour-coded according to the sign of Q ra d (dark blue and 
yellow correspond to <3 ra d < and Q ra d > 0, respectively). Lower panel: Trajectories are colour- 
coded according to the value of T — Tre, the deviation of the gas temperature from the radiative 
equilibrium temperature. In both panels, the fluid elements emanate from the left-hand-side and 
travel towards the right. 



blue means Q ra d < 0, yellow means Q ra d > 0). The motion of the fluid element before it emerged 
at the surface is essentially adiabatic, with s = 6.1R*, where R* is the universal gas constant. As 
it approaches z « 0, radiative cooling leads to a steep drop of s and T with height. Accompanying 
this cooling is a simultaneous compression of the fluid element, so that it becomes temporarily 
denser than the fluid beneath it. Thereafter, the fluid clement continues to ascend and the density 
decreases again. When it has reached z = 0.1 Mm, the fluid element is cooler than Tre and 
it begins to be heated radiatively. This radiative heating leads to a rise in entropy, while the 
temperature continues to decrease because the effect of expansion overwhelms that of radiative 
heating. Finally, at z = 0.28 Mm, the fluid element overturns and begins its descent. At this 
moment, the fluid element has a gas temperature below Tre, so that it continues to be heated 
radiatively even though its temperature rises by compression. Only at z = 0.26 Mm (w 20 km 
below the apex of the fluid element's trajectory) is the temperature of the fluid element above Tre 
and Q ra d becomes negative. 

The profiles of s, T and g in Fig. [8] exhibit hysteresis. There are two reasons for this behaviour. 
Firstly, there is a lag between when the fluid clement overturns and when Qrad switches sign (see 
previous paragraph). Secondly, the asymmetry in the speed of upflows (slower) and downflows 
(faster) means that the cumulative heating experienced by the fluid element over its journey up- 
wards in the photosphere is greater than the accumulative radiative cooling it undergoes on the 
return journey. These two factors contribute to the temperature and entropy enhancements of the 
intergranular boundaries (downflows) with respect to the granules (upflows). 

In order to test whether dissipative effects (viscous and thermal, e.g. in shocks) play a role in 
causing the temperature enhancement in the intergranular lanes in the upper photosphere, we have 
reconstructed the profiles of s along the trajectory neglecting hyperdiffusivity effects. By keeping 
track of Q r ad(i), g(t) and T(t) along the trajectory of a fluid element, we can reconstruct the time 
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Figure 8: Profiles of specific entropy s, temperature T and density g along one of the fluid element 
trajectories shown in Fig. As the fluid element initially emerges at the base of the photosphere 
(z « 0), it has relatively high entropy, high temperature and low density. The colour-coding 
indicates the sign of Q ra d experienced by the fluid element along its trajectory (dark blue means 
Qrad < 0, yellow means Q ra d > 0). Above z — 140 km, the fluid element is hotter during descent 
than during its ascent. This is a result of radiative heating of the fluid element in the photosphere. 
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Figure 9: Reconstruction of specific entropy profile. The black solid curve shows the entropy profile 
as obtained by following the fluid clement in the simulation (same as in the upper panel of Fig. [8| . 
The red dashed curve shows the reconstructed entropy profile, s rccon , calculated by integrating Eq. 

history of s for a fluid element. Given these profiles and the initial entropy sq of the fluid element, 
Eq. ([7]) can be integrated to yield the reconstructed entropy profile s TCCOn (t). This quantity is 
plotted in Fig. 03 together with the entropy profile shown in Fig. [HI The good match of the two 
profiles in the layers z > tells us that radiation is indeed the cause of hysteresis in the entropy 
and temperature profiles. The reconstructed profile s rcC on begins to deviate from the tracer profile 
once the fluid element approaches the narrow turbulent downflow jets in the convection zone. Here 
the radiative effects diminish and dissipative effects described by the hyperdiffusive terms become 
important. 

4 Conclusion 

The reversed granulation is a reversal of the pattern of intensity contrast between intensity maps 
in the continuum and intensity maps in spectral lines forming in the middle to upper photosphere. 
It has long been established that this reversal of the intensity fluctuations is a result of the sign 
change of the horizontal temperature fluctuations (see Fig. [5]). Analysis of our 3D numerical 
model of the near-surface convection zone and photosphere confirms the existence of this reversal 
in the temperature fluctuations above a height of z = 130 — 140 km, in agreement with previous 
observational findings. 

By following in detail the time histories of fluid elements overshooting into the photosphere, we 
were able to identify the processes that determine the horizontal temperature fluctuations in the 
photosphere. We find that changes in the temperature of fluid elements in the photosphere result 
from adiabatic expansion and radiative heating. Whereas the former decreases the gas temperature 
of a fluid element, the latter increases it towards the radiative equilibrium temperature Tre- The 
ratio of the associated characteristic timescales tells us that, for a fluid element rising in the 
optically thin layers, the decrease of temperature due to expansion is only partially cancelled by 
radiative heating. This disparity allows a fluid element to depart from radiative equilibrium as 
it penetrates into the photosphere. Due to this deviation from Tre, a fluid element overturning 
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in the optically thin layers still experiences some radiative heating as it begins to descend. This, 
together with the fact that downflows have higher vertical speed than upflows, leads to a hysteretic 
behaviour of the thermodynamic history of the fluid element (see Fig. [SJ. 

In various guises in previous works, the following explanation has often been invoked to explain 
the origin of the reversal of temperature fluctuations: since the photosphere is subadiabatically 
stratified, material rising from granules into this stable layer has lower temperature that its 'mean 
surroundings'. This view is problematic in two respects. Firstly, it implicitly assumes that ma- 
terial penetrating into the photosphere does so in an adiabatic fashion. Secondly, it is not at all 
clear what are the 'mean surroundings' of a fluid element embedded in convective flow consisting 
of disjointed upwellings separated by an intergranular network. Our present efforts highlight the 
importance of radiative energy exchange in the dynamics of the photosphere and helps dispel pos- 
sible misconceptions regarding the origin of the reversed granulation in the photosphere. 
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